If sterile neutrinos exist, how can one determine the total solar neutrino fluxes? 
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The 8 B solar neutrino flux inferred from a global analysis of solar neutrino experiments is within 
11% (lcr) of the predicted standard solar model value if only active neutrinos exist, but could be 
as large as 1.7 times the standard prediction if sterile neutrinos exist. We show that the total 8 B 
neutrino flux (active plus sterile neutrinos) can be determined experimentally to about 10% (lcr) by 
combining charged current measurements made with the KamLAND reactor experiment and with 
the SNO CC solar neutrino experiment, provided the LMA neutrino oscillation solution is correct 
and the simulated performance of KamLAND is valid. Including also SNO NC data, the sterile 
component of the 8 B neutrino flux can be measured by this method to an accuracy of about 12% (lcr) 
of the standard solar model flux. Combining Super-Kamiokande and KamLAND measurements and 
assuming the oscillations occur only among active neutrinos, the 8 B neutrino flux can be measured 
to 6% (lcr); the total flux can be measured to an accuracy of about 9%. The total 7 Be solar neutrino 
flux can be determined to an accuracy of about 28% (lcr) by combining measurements made with the 
KamLAND, SNO, and gallium neutrino experiments. One can determine the total 7 Be neutrino flux 
to a la accuracy of about 11% or better by comparing data from the KamLAND experiment and the 
BOREXINO solar neutrino experiment provided both detectors work as expected. The pp neutrino 
flux can be determined to about 15% using data from the gallium, KamLAND, BOREXINO, and 
SNO experiments. 



I. INTRODUCTION 



We describe in this paper analysis procedures that can answer two of the most important questions of neutrino 
research. How can one determine the total solar neutrino fluxes ( 8 B, 7 Be, and pp) for comparison with solar model 
predictions? How can one determine the sterile contribution to the total solar neutrino fluxes? Our answers allow 
for the possibility of an arbitrary mixture in solar neutrino oscillations of active and sterile neutrinos, but require 
the correctness of the LMA solution of the solar neutrino problems and careful attention to all the sources of error 
(theoretical as well as experimental) [JJ 

We focus first on determining total fluxes by comparing charged current (CC) observables measured in different 
experiments; this method yields results as independent as possible of uncertainties due to the presence of sterile 
neutrinos. We then describe how similar techniques can be applied to determine total solar neutrino fluxes using a 
CC experiment plus a neutrino-electron scattering experiment [or a neutral current (NC) measurement], which yields 
results that depend more on the sterile neutrino mix but which can nevertheless be relatively accurate. 
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1 This paper was originally written and posted on the electronic archive (hep-ph) before the announcements of the recent 
SNO results |lj^| and the improved SAGE measurement of the gallium rate || and also before our paper was submitted for 
publication. We have included in the analysis reported in this version of the paper, which we are submitting for publication, the 
recent SNO and SAGE measurements. The ideas with respect to the 7 Be and 8 B neutrinos are unchanged and the numerical 
results have not been affected significantly, but the present version is more up-to-date with respect to the input data. We have 
also added, inspired by the SAGE discussion, a detailed analysis of what one can learn about pp neutrinos before there is a 
dedicated experiment to measure just the pp neutrino flux. 
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The numerical values we estimate for the expected precision with which different quantities can be measured 
rely upon simulations of the performance of the relevant experiments. Therefore the accuracies that we quote are 
illustrative; the actual accuracies that are obtainable can only be determined once the experimental uncertainties are 
known. 

A. Flavor changes occur 

Neutrinos change flavors as they travel to the Earth from the center of the Sun. This flavor change was seen directly 
by the comparison of the Sudbury Neutrino Observatory (SNO) measurement [Q of the charged current reaction for 
8 B solar neutrinos with the Super-Kamiokande measurement Q of the neutrino-electron scattering rate (charged 
plus neutral current). Even more clearly, flavor change has been demonstrated by comparing the neutral current 
measurement by SNO with the SNO CC measurement [|J. The conclusion that flavor changes occur among solar 
neutrinos, if based solely upon the comparison of the SNO and Super-Kamiokande event rates, is valid statistically 
at about the 3.2cr confidence level §]• The neutral current measurement of SNO increases the significance level 
for flavor changes among solar neutrinos to the 5.3<r confidence level. 

The SNO and Super-Kamiokande results demonstrate simply that new physics is required to resolve the long- 
standing solar neutrino problem pJOj ] , i.e., to understand the origin of the discrepancy between the predictions of the 
standard solar model |fl[] and the observed solar neutrino event rates J|J|,[T^-[h| . If one includes the results of the 
Chlorine @, Kamiokande ||, SAGE @, GALLEX JT|, and GNO ||] experiments together with the SNO (CC) 
and Super-Kamiokande results, then the combined measurements require j^] new physics at 4.0a and, if the relative 
temperature scaling of the 7 Be and 8 B neutrino production reactions is taken into account, at 7 ' Aa. Helioseismological 
measurements confirm the predicted sound speeds of the Standard Solar Model to better than 0.1% and show that 
stellar physics cannot account for the discrepancies between standard predictions and the observed solar neutrino 
rates 0. 

B. Current knowledge of the 8 B solar neutrino flux if only active neutrinos exist 

The combination of the charged current (CC) and the charged plus neutral-current measurement with Super- 
Kamiokande has been used by several groups to determine the flux of active 8 B neutrinos independent of the 
solar model. These model-independent determinations of the active flux exploit the similarity between the response 
functions in the SNO and Super-Kamiokande detectors 18 In addition, if one includes all the experimental 



data in a global oscillation solution in which the 8 B flux is a free parameter, one obtains a similar (but slightly smaller) 
allowed range for the 8 B neutrino flux pp|-p2]|. The measurement of the NC rate by SNO [Q provides an independent 
determination of the active 8 B neutrino flux. 

All of the analyses yield the same result: if electron neutrinos oscillate into only active neutrinos, then the total 8 B 
neutrino flux is in excellent agreement with the flux predicted by the Standard Solar Model. 

This close agreement of the active 8 B neutrino flux with the total flux predicted by the Standard Solar Model 
(SSM) |ll],|22| is, if the flux of sterile neutrinos is small, an important confirmation of the quantitative theory of 
stellar evolution. We summarize below the current best-estimates and the associated la uncertainties for the active 
8 B neutrino flux, </> ac tivo( 8 B). 

• Standard solar model (BP00) prediction Q: Q 

0activc( 8 B) = 5.05 x lOWs-^l+g;*?). (I) 

• Comparison of SNO and Super-Kamiokande event rates (cf . esp. Ref . Q ) : 

</>activc( 8 B) = 5.44 x 10 6 cm 2 s^(l ± 0.18). (2) 



2 The recently measured low-energy cross section factor reported by Junghans et al. p3] is currently being reinvestigated. 
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SNO neutral current measurement (assuming undistorted spectrum) Q . 

0activc( 8 B) = 5.09 x 10 6 cm 2 s _1 (l ± 0.12). (3) 

Global neutrino oscillation analysis (see Table this paper and Ref. E2] ) : 



^active 



( 8 B) = 5.40 x 10 8 cm 2 s -1 (l ± 0.075). (4) 



The agreement, summarized in Eqs. (jl|)-(^), between the SSM calculated flux and the measured active flux suggests 
that the sterile neutrino contribution to the 8 B neutrino flux may be small. In this paper, we ignore this tempting 
suggestion and instead concentrate on developing methods to determine experimentally the total 8 B and 7 Be neutrino 
fluxes emitted by the Sun, independent of the active-sterile mixture (for a discussion of earlier investigations of sterile 
neutrinos see Refs. 24,2j|). 

If we want to understand the particle physics implications of solar neutrino research, we must determine if sterile 
neutrinos are present in the solar neutrino flux . Moreover, the original-and still valid-goal of solar neutrino research 
was Jul to compare solar model predicted and experimentally measured (total) solar neutrino fluxes. 



C. What can one do if sterile neutrinos exist? 



What is the situation if sterile neutrinos exist? The total flux of 8 B neutrinos could in this case be much larger than 
the standard solar model prediction; a major fraction of the total flux that reaches the Earth could arrive in a form 
that is not detected in solar neutrino experiments. The existing data disfavor (at 5.4cr CL) oscillation into purely 
sterile neutrinos. Nevertheless, a large sterile component is allowed J2^,|26| if oscillations occur into a combination of 
active and sterile neutrino states (see Refs. |^],^8| for a description of the formalism adopted here). The flux of sterile 
neutrinos could in principle be large enough to destroy the apparently excellent agreement between the flux predicted 
by the SSM and the true flux of 8 B neutrinos [which is assumed to be pure active neutrinos in the comparison shown 
above in Eqs. (|) -(§)]. 

A measurement of the total 8 B solar neutrino flux, including the sterile component (if any), will provide information 
that is important for astrophysics and for particle physics. The motivation for investigating sterile neutrinos is not 
dependent upon the LSND |^9| results that might suggest the existence of sterile neutrinos. Of course, the LSND 
results are not supported by other experimental results (cf. Ref. [j30| ) nor by theoretical predictions that have been 
precisely confirmed (like the helioseismological verifications of the standard solar model) , as is the case for the inference 
of flavor change based upon the SNO-Super-Kamiokande comparison. 

Fortunately, the KamLAND reactor neutrino experiment [3l| ] , when combined with the SNO measurement of the 
CC flux, is capable of providing a precise determination of the total, i.e. the active plus the sterile, 8 B neutrino flux. 
We assume throughout this paper the correctness of the currently favored large mixing angle (LMA) solution to the 
solar neutrino problem. If the LMA solution is not correct, then all global analyses of the available solar and reactor 
data indicate that either the mass difference, Am 2 , or the vacuum mixing angle, 9, will be too small to produce a 
measurable effect in the KamLAND experiment (see e.g. Refs. [p|— p| j2C| , |2l] , |3l| , p2| | ) . In this case, KamLAND will be 
unable to provide the information required to determine the total 8 B neutrino flux. 

Here is the basic physical idea of the method we propose for measuring the total 8 B solar neutrino flux. For the 
KamLAND J3lj experiment, one will know accurately the flux of anti-neutrinos from the 17 reactors that contribute 
significantly to the measured anti- neutrino events. From measurements of the total event rate and the energy spectrum 
induced by the surviving z7 e , the KamLAND experimentalists can determine with precision (3l],|3^-|35| the anti- neutrino 
propagation parameters, Am 2 and tan 2 9. Both the KamLAND measurement and the CC SNO measurements are 
disappearance experiments for neutrinos (or anti-neutrinos) of similar energies. For the CC measurement made with 
SNO, one does not know the total 8 _B neutrino flux created in the Sun. But, assuming conservation of CPT, one can 
use the propagation parameters Am 2 and tan 2 9 determined by KamLAND and the measured (by SNO) CC rate to 
solve for the flux that gives the observed result. Summarizing, for the KamLAND experiment one knows the total 
flux but not the propagation parameters, which are measured. For the SNO CC experiment, one will know (from 
KamLAND) the propagation parameters and therefore can measure the total flux^|. 



3 The method described here is of course more general than the specific application to the KamLAND and SNO experiments. 
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FIG. 1. Solar neutrino allowed region compared with simulated KamLAND allowed region. The figure shows the currently 
allowed regions of the solar neutrino oscillation parameters; the contours of equal CL are labeled at la, 2a and 3a. This global 
solution was obtained assuming pure active neutrino oscillations and using all the available solar and reactor data. We include 
the recent SNO results The rates from the GALLEX/GNO and SAGE ||Q experiments have been averaged to 

provide a unique data point (72.4 ±4.7 SNU). Some technical improvements regarding neutrino cross sections and correlations 
of errors were included in the analysis (see the Appendix). The two much smaller allowed regions, placed symmetrically with 
respect to the line at tan 2 9=1, represent the allowed regions, at la, 2a, and 3a, that are obtained from a simulation of 
what may be achievable with the KamLAND reactor experiment. The best-fit point for the KamLAND simulation is assumed 
to be the same as the best-fit point for the global solar neutrino oscillation solution, namely, purely active neutrinos with: 



Am 2 



5.0 x 10~ b eV 2 , tan 2 9 = 0.42. 



Figure hi shows the results of a refined global solution for the solar neutrino oscillation parameters that was made 
(see Sec. O and the Appendix) using all the available solar and reactor data. The figure displays the allowed solar 
neutrino oscillation contours at la, 2a, and 3<r. The results are obtained by the procedures described most recently 
in Ref . [£2| , where we have used in the present paper the analysis strategy (a) (of Ref. ) including the 1496 day 
Super-Kamiokande data sample |37j] as well as the SNO CC, NC and day- night observations 0^]. We have also made 
some improvements (see the Appendix) in the treatment of the uncertainties in the neutrino cross sections and in the 
correlation of errors. 

For comparison, we also show in Fig. |j] the small size of the expected allowed region for KamLAND if this reactor 
anti- neutrino experiment observes a signal corresponding to the current best-fit point of the solar neutrino analysis. 
In calculating the KamLAND allowed region, we have made a conservative estimate (see Sec. IIIC), following the 
principles discussed in Refs. ||^,|3£j. Figure [j] shows clearly that KamLAND has the potential for making a precise 
measurement of the solar neutrino oscillation parameters, provided that the LMA is the correct oscillation solution. 

The only complication involved in determining the total 8 B flux from a comparison of the KamLAND and SNO CC 
measurements results from the fact that for the favored Large Mixing Angle (LMA) solar neutrino oscillation solution 
matter effects in the Sun and the Earth can be significant. Matter effects are unimportant for the KamLAND reactor 
experiment. The role of matter effects in solar neutrino experiments depends somewhat upon the a priori unknown 
active-sterile mixture, which introduces a calculable uncertainty in the inferred total 8 B neutrino flux. 

In summary, the combination of the SNO CC neutrino measurement and the KamLAND anti-neutrino measurement 
will determine the total (active plus sterile) flux, 0totai( 8 B), of 8 B solar neutrinos. By subtracting the previously 



In order to determine the total flux, it is sufficient that one measure a set of observables that do not depend upon the solar 
neutrino flux (in this paper, the measured quantities in the KamLAND experiment) and a quantity that does depend upon the 
solar neutrino flux (here, the CC rate in SNO). 
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determined flux (see above), </> ac tive( 8 B), from </>totai( 8 B), one can determine the flux of sterile solar neutrinos. 

Using similar reasoning, we shall also show that the combined Super-Kamiokande and KamLAND measurements 
can be analyzed to yield an accurate value for the total 8 B flux, although in this case the results are somewhat more 
sensitive to the active-sterile admixture. 



D. Be solar neutrinos 



The flux of 7 Be solar neutrinos, which in the SSM is predicted Q to be 0( 7 Be act ivc) = 4.77 x 10 9 cm 2 s _1 (l±0.10), 
can be determined in a model independent way from a global analysis of the solar neutrino data assuming 
only active neutrino oscillations. For example, the latest analysis by Garzelli and Giunti Q yields 0.02 < 
(/>( 7 Be)M 7 Be,SSM) < 1.15 at 99% CL. 

We shall also show in this paper that one can extract the value of the 7 Be neutrino flux from measurements of 
the gallium solar neutrino experiments, GALLEX, SAGE, and GNO and the results of the SNO and KamLAND 
measurements. The value of the 7 Be flux that will be derived in this way is relatively insensitive to the assumed 
neutrino oscillation parameters, although it does depend somewhat on the assumed contributions of the CNO, pep, 
and pp neutrino fluxes which we adopt from the standard solar model. The constraint provided by the Chlorine 
experiment is not very useful for determining the total 7 Be neutrino flux. 

One can obtain an independent measurement of the total Be solar neutrino flux by comparing data from the 
KamLAND experiment with data from the BOREXINO [^(J solar neutrino experiment. If both the KamLAND and 
the BOREXINO detectors work as expected, then this method will be more accurate than the methods involving the 
gallium and chlorine radiochemical detectors. 



E. Appendix: just for aficionados 



The determination of the total solar neutrino fluxes, and even more so the determination of the sterile components 
of these neutrino fluxes, requires precision in both the experimental measurements and the theoretical calculations 
and analyses. We present in the Appendix a refined discussion of the theoretical errors, and their correlations, for the 
absorption cross for the gallium and Chlorine solar neutrino experiments. 



F. Outline and suggested reading strategy 



The outline of this paper is as follows. In Sec. [n], we describe the current experimental knowledge of the 8 B solar 
neutrino flux. Our results are summarized in Table Q both for the special case of oscillations between purely active 
neutrinos and for the general case of oscillations between electron neutrinos and an active-sterile neutrino admixture. 
We limit our analysis to the allowed LMA region of solar neutrino oscillations. We show in Sec. Ill how one can 
use the CC measurements with SNO and KamLAND to determine an accurate total 8 B solar neutrino flux including 
experimental and theoretical uncertainties and the possibility of an appreciable active-sterile admixture. We switch 
to the 7 Be flux in Sec. |y| and evaluate how well one can determine the total 7 Be solar neutrino flux by also using the 
results of the gallium (GALLEX, SAGE, and GNO) solar neutrino experiments or the Chlorine experiment. In Sec. [v|, 
we investigate how well the total 8 B and 7 Be neutrino fluxes can be determined using the combined measurements of 
KamLAND and v — e scattering observed in the Super-Kamiokande (Sec. |V A| ) and BOREXINO (Sec. |VB| ) detectors. 
We show that even in the presence of active-sterile admixtures the total 'Be solar neutrino flux may be measured 
with relatively high accuracy by comparing results from the KamLAND and the BOREXINO experiments. We 
describe and analyze in Sec. VI three strategies for determining the total pp solar neutrino flux in the absence of a 
dedicat ed e xperiment that measures separately the pp neutrinos. We summarize and discuss our principal conclusions 
in Sec. |yj| 

We urge the reader to turn first to Sec. VII and read there the summary and discussion of our main results and 
their implications. The rest of the paper can then be understood more easily. 



II. PRESENT KNOWLEDGE OF THE S B NEUTRINO FLUX 



We generalize in this section the determination of the 8 B neutrino flux to the case in which v e oscillates into a state 
that is a linear combination of active (v a ) and sterile (y s ) neutrino states, 
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v e — > cos n]v x + sin 77 v s , (5) 

where 77 is the parameter that describes the active-sterile admixture. This admixture arises naturally in the framework 
of 4-j/ mixing pl[ ). The total 8 B neutrino flux can be written 

<M 8 B) tot al = <t>{Ve) + <t>{Vx) + <f>{Vs), (6) 

where 4>{v s ) — tan 2 77 x 4>{v x ). Clearly, the larger the sterile component, the larger the value of </>( 8 B) to tai that is 
inferred from the experimental data. 

We have performed a global analysis of the solar neutrino data treating the total 8 B neutrino flux as a free parameter. 
The details of the analysis procedure are the same as those used in Ref. p2j except where we explicitly state otherwise. 
We concentrate here on the LMA region, 0.1 < tan 2 9 < 10, 10~ 5 < ArrP/eV 2 < 10 -3 . 

To take account of the possibility of oscillations into sterile neutrinos, we determine the allowed regions in the 
parameter space defined by Am 2 , tan 2 9, and a third parameter, cos 2 r), that is defined by Eq. (^|). It is convenient to 
introduce the dimensionless parameter, 

„ _ 0( 8 B)total 

lB 0( 8 B) SSM ' {,) 

where /b is the total 8 B neutrino flux in units of the predicted standard solar model flux. The allowed range of /b in 
the three dimensional space of neutrino parameters, Am 2 , tan 2 9, and cos 2 77, is determined by the equation 

X 2 (f B ) < xLn + A X 2 (1,CL). (8) 

Here Ax 2 (l, CL) is the change in \ 2 that corresponds to a specified confidence limit (CL) for one degree of freedom. 
The computed values of \ 2 are minimized for each value of /b with respect to Am 2 , tan 2 9, and cos 2 77. 

Table [I] shows the currently allowed range for fg for both the more general case where a mixture of active and 
sterile neutrinos is assumed and for the more conventional case in which only active neutrinos are considered. For 
purely active neutrinos, the la range is 

/ B , active = 1.07 ±0.08, (9) 
and for an arbitrary mixture of active and sterile neutrinos, the la range is 

/b, active — sterile = 1.07_ q Q8 ' (10) 

The result shown earlier in Eq. for purely active neutrinos is taken from Table ffl. 

How does the possible existence of sterile neutrinos affect the allowed range of 8 B neutrino fluxes? We can calculate 
the dependence of the allowed range of Jb upon cos 2 77 with the aid of the inequality 

X 2 (/ B ,cos 2 7 7 ) < X 2 lin + A X 2 (2,CL). (11) 

Here A% 2 (2, CL) is the change in x 2 f° r a specified CL that corresponds to two degrees of freedom ; the computed 
values of \ 2 are minimized at each point with respect to Am 2 , and tan 2 9. 

Figure U shows the range of /b as a function of the active-sterile admixture, cos 2 77, that is obtained from Eq. (|ll|) for 
the la, 2a and 3a allowed regions. The allowed regions are defined respect to the global minimum, which corresponds 
to purely active oscillations with Am 2 = 5.0 x 10~ 5 eV 2 , tan 2 9 — 0.42, and f B = 1.07. 

Although pure sterile oscillations are forbidden at the 3 a CL (cf. Fig. ||), a large sterile admixture in the solar 
oscillations is still allowed. In fact, with the currently available data, the largest allowed value at 3c of the sterile 8 B 
neutrino flux corresponds to /b = 2.2 and cos 2 77 = 0.3 (for 2 d.o.f.). For this extreme case, we find 

TABLE I. The allowed range of the total 8 B neutrino flux. The table presents the allowed range of }b [defined by Eq. (Q)] 
that was found in a global solution of all the currently available solar and reactor neutrino data. The second column gives the 
allowed range of /b for an arbitrary mixture of active and sterile neutrinos and the third column gives the range for the case 
where only active neutrinos are considered. The results shown were obtained using Eq. (pi). 
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FIG. 2. The dependence of the inferred 8 B flux on the active-sterile admixture. The figure shows, as a function of the 
active-sterile admixture, i.e., cos 2 i], the allowed range of the 8 B solar neutrino flux at la, 2a and 3a CL. The star indicates 
the global best fit point for all the currently available solar and reactor data; the star lies at /s = 1.07 and r\ = 0.0 (purely 
active neutrinos). 



/] 



B. sterile max 



l.i, 



(12) 



The quantity /b, sterile max that appears in Eq. (|l2|) is defined, analogous to Jb in Eq. (0), by the relation 
/b, sterile max = 0( 8 B) s tcriie max/0( 8 B) S sM- The maximum value of / B , sterile max is as large as the sum of the ac- 
tive 8 B neutrino fluxes (f x + f e = 1.1, where f e — 0.348, see Ref. for this special caseQ. 

What is the maximum allowed sterile contamination of the 8 B solar neutrino flux? Minimizing x 2 for the global 
solution with respect to Am 2 , tan 2 9, and we find that the allowed range of cos 2 rj satisfies 



0.75(0.40) < cos 2 77 < 1.0 



(13) 



at ler(or 3a). 



III. HOW CAN WE DETERMINE THE TOTAL S B NEUTRINO FLUX USING CC REACTIONS? 

In this section, we will show how one can determine the allowed range of the total 8 B neutrino flux using the results 
of the KamLAND reactor neutrino experiment and the SNO CC solar neutrino experiment. We shall also estimate 
the accuracy with which one can determine the total 8 B flux. 

Since we consider here only CC reactions that result from disappearance experiments, the only difference between 
the role of active neutrinos and v T and sterile neutrinos ^sterile arises from matter effects in the Earth and in the 
Sun. Since sterile neutrinos do not interact with matter, the effective potential for the v e —v a evolution in matter is 
V cs = V c — V s = Vcc + ^nc, since V s = 0. The effective potential V es is approximately half the potential for v e -v a , 
V ca = V c — V a — Vcc 1 where V a is the potential for the active neutrinos and v r . (The difference is exactly half for 
a medium with equal number of neutrons, protons and electrons because V a = Vmc = —GFN n /y2 ~ — Vbc/2 with 
Vcc — V^GfiVe.) We shall evaluate the expected dependence of the inferred total 8 B flux on the admixture of sterile 
neutrinos [cf. Eq. (H)]. 



4 For 2 d.o.f, the maximum allowed value of fs is 2.2 at 3<r, but is 1.7 for 1 d.o.f, see Table |. We have given in the Abstract 
the maximum value for 1 d.o.f. 



7 



In Sec. [II A we present the for mulae that are used to determine the 8 B flux with the aid of the KamLAND and 
SNO CC experiments and in Sec. Ill B w e illus trate the effect on the inferred 8 B flux of the maximum allowed (at 
la) sterile admixture. We estimate in Sec. Ill C the precision with which the 8 B flux can be determined including all 
the principal known sources of uncertainties. 

The reader who is interested in how well we can determine the 8 B flux, but does not need to know the details of 
the proce dure , can get the main results by glancing at Fig. || and Fig. || and Table ||. 

In Sec. V A, we investigate how well the total 8 B flux can be determined using KamLAND in combination with the 
v — e scattering experiment, Super-Kamiokande. 



A. Relations that determine the total flux 



Suppose KamLAND observes a signal that corresponds to LMA T> e oscillations with parameters (Am 2 , tan 2 ). We 
assume the validity of the CPT theorem so that constraints on anti-neutrino oscillation parameters obtained from 
the KamLAND experiment apply to solar neutrino experiments. We can then extract the 8 B neutrino flux from the 
following relation: 



7-)CC.cxp -. 
^SNO v £_ 



i?||M (P ee (Am 2 ,tan 2 0)) S No' 



(14) 



where 



#1no = / dE u SSM ( 8 B, E u )a e {E v ) = 2.87i?^ Q cxp (15) 

is the CC rate for the SNO experiment Q that is predicted 0,^2| by the standard solar model in the absence of 
oscillations and (P ee (Am 2 , tan 2 0))sno is the average survival probability for electron-flavor neutrinos created in the 
Sun. Also, E v is the neutrino energy and cr e is the weighted average v e -d interaction cross-section, including the 
experimental energy resolution function, Res(T, T"), where T(T") is the measured (true) recoil kinetic energy of the 
electron. Thus 

/•T max r T mMC '(E u ) j r-p rp,\ 

a e (E u ) = J t dT jf dT' Res(T, T') ^g,' j . (16) 

The lower limit, T m j n , in the integral in Eq. ( |l6| ) is taken here to be the threshold used by the SNO Collaboration in 
Ref. (T m ; n = 5 MeV). The calculated value for the CC rate is not sensitive to the assumed value of T max , as long 
asT max > 17MeV. 

The energy-averaged survival probability, (P ee (Am 2 , tan 2 #))sno for v e at SNO can be computed using the propa- 
gation parameters, (Am 2 , tan 2 0) , observed at KamLAND. Thus 



JdE v SSM ( 8 B, E v )oAE v )P ee {E v , Am 2 , tan 2 6) 



{P ee {Am ,tan 8)) S no = ~ ^§sm • ( 17 ) 



L SNO 



B. Illustrative dependence of total flux upon active-sterile admixture 

How much does a sterile neutrino admixture affect the inferred total 8 B neutrino flux? The dominant dependence 
on the sterile admixture arises from matter effects within the Sun for larger Am 2 and within the Earth for smaller 
Am 2 . 

Figure || shows the isocontours of /b in the LMA region for the pure active case (thicker lines) and for the active- 
sterile case cos 2 r\ = 0.75 (thinner lines). The isocontours are not very different for the two cases. Moreover, the 
dependence on Earth matter effects can be avoided experimentally by using only the daytime CC measurement, once 
sufficient statistics are available. The robustness of the inferred 8 B solar neutrino flux can be tested by comparing 
the 8 B flux inferred using only daytime CC measurements with the flux that is inferred when nighttime data (with 
corrections for Earth matter effects) are added to the daytime data. 

In Fig. ||,we also show the 3a LMA contour (the dotted contour) obtained by the global analysis of the solar neutrino 
data for the purely active case . Within the 3cr LMA region, the maximum difference between the value of /b inferred 
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FIG. 3. Isocontours for the total B neutrino flux. The figure compares isocontours for the B flux assuming purely active 
neutrino oscillations (thicker lines) with the flux that would be inferred for a 75% active-25% sterile admixture (thinner lines). 
The results refer to a hypothetical comparison of measurements from the KamLAND reactor experiment and the SNO CC 
experiment. We also show (dotted contour) the 3<r allowed region obtained by a global fit to all of the allowed solar and reactor 
data (cf. Fig. [j]). 



allowing for possible sterile neutrinos and the value obtained assuming only active neutrino oscillations is +0.9% and 
—3.5%. The dependence upon the active-sterile mixture could become negligible if the correct Am 2 lies in the lower 
part of the LMA region and only daytime data is used from the SNO CC measurements. 

We conclude from Fig. [| that existence of sterile neutrinos will not prevent an accurate measurement of the total 
8 B neutrino flux. 



C. How accurately can the total B flux be determined? 

What is the overall precision expected in the determination of the total 8 B flux? From Eq. (|l4|) we can derive the 
anticipated precision as 



*(/b) 



<t{FL 



CC.exp 
SNO 



R 



CC.cxp 
SNO 



a(R 



SSM\ 
SNO J 



pSSM 
-"-SNO 



ee/SNO, KamLAND 



(-Pee) SNO, KamLAND) 



= (5 (/B)s N o.cxp + <*(/b) 



/SNO, Cross Section 



<$(/] 



BJsNO, KamLAND' 



(18) 



where the la errors are combined quadratically. 

The current value of the first term in Eq. (jig), the la uncertainty (statistical and systematic) of the measured CC 
rate in SNO, is Jl| ^(/b)sno cxp = 6.15%. This uncertainty will undoubtedly decrease as the results of analyzing more 
SNO CC data are reported. 

The second term in Eq. (18) represents the uncertainty in the u e - 2 H absorption cross section. Much progress has 
been made recently in evaluating this cross section, see e.g. Refs. fFl| , [42| ], which has led to an estimate ~ 1% for 
the cross section uncertainties other than radiative corrections. No definitive calculation has yet been made of the 
radiative correction for the CC reaction, but a reasonable estimate |yj is that the cross sections given in Ref. plj | 
might be increased by 2%. We adopt here a conservative uncertainty of <5(/b)sno,c.s. = 2% ; the precise value chosen 
for <5(/b)sno.c.s. is not very important at this stage since other uncertainties are dominant. The hep neutrino flux 
contribution to A/b is negligible for our purposes j||5||ll[] 

A detailed simulation is required to estimate cr((P ee )sNO, KamLAND), i.e., the uncertainty in the average electron 
neutrino survival probability for 8 B solar neutrinos observed in the SNO CC experiment as will be determined by 
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TABLE II. Values of /b and associated uncertainties obtainable from the SNO CC and KamLAND experiments. The table 
presents the best-fit values for /b (the total 8 B neutrino flux divided by the predicted standard solar model 8 B neutrino flux) 
and associated uncertainties for a representative set of possible oscillation parameters. We have used 5(/b)sno,cx P = 6.15% 
and <5(/b)sno,cs = 2%. We consider active-sterile neutrino admixtures permitted by the currently allowed global oscillation 
solution [see Eq. (|l2|)]. 
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+6.9 
-8.1 
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future KamLAND measurements. Here is how we estimate this uncertainty. We generate the expected KamLAND 
signal for a fine grid of points, (Am 2 , tan 2 9) that spans the space of the allowed oscillation parameters determined 
from solar neutrino experiments. For each grid point, we obtain the allowed region by a % 2 analysis. We use statistical 
errors corresponding to three years of data taking at KamLAND, observing anti- neutrinos from reactors working at a 
constant 78% of the maximal power. To be conservative, we also assumed a neutrino energy threshold of 3.5 MeV, in 
order to ensure that the effects of natural radioactivity would be small. More details on the KamLAND experiment 
can be found in Ref. |Q; details regarding the neutrino cross sections, statistical procedures, and reactor fluxes used 
in the present paper are described in Refs. ]3^j39[ |. 

In computing the inferred values of /b, we take account of the fact that there could be a significant component of 
sterile neutrinos in the incident 8 B solar neutrino flux. We therefore consider all active-sterile admixtures permitted 
by the global oscillation solution shown in Fig. [l]. The numerical constraint on the currently allowed admixture is 
given in Eq. ([i~3"|). 

In principle, for each simulated point there are two allowed KamLAND regions, one around Am 2 and tan 2 9 and 
another around Am 2 and tan 2 ( : | — 9). We discuss here only the range of parameters within the first octant for 
the mixing angle since global solar neutrino solutions show Jf^-|8|,p0|- p2| , ^2[ that the first octant is preferred. For 
each simulated KamLAND allowed region centered on a specific (Am r\ tan 2 9), we compute cr((P ee )KamLAND)- We 
repeated this procedure for a grid of 81 x 41 points in the range 0.1 < tan 2 9 < 10, 10~ 5 < Am 2 < 1 x 10~ 3 . The 
estimated uncertainty in (-Pee)sNO, KamLAND varies with the grid point (Am 2 , tan 2 9). 

Table [0| presents the best-fit value of /b, the uncertainty in cr((P ee )sNO, KamLAND), and the total expected uncer- 
tainty in inferring Jb from the combined SNO CC and the KamLAND neutrino reactor measurements for a represen- 
tative set of possible results for Am 2 and tan 2 9. In all the cases shown in Table ||, we have used £>(/b)sno,cx P = 6.15% 
and <5(/B)sno,cs = 2% (see previous discussion). The uncertainty contributed by the possibility that sterile neutrinos 
exist, i.e., r\ ^ 0, is rather modest; cr((P ee )sNO, KamLAND) is typically reduced by ~ 1% from the value shown in 
Table [l| 

Figure ^ shows the 9% and the 14% contours for the maximum percentage deviation from the best-fit /b value. 
To provide a context, the figure also displays the ler, 2cr and 3<r allowed LMA regions obtained by a global fit to the 
available solar and reactor data. Within almost all of the current 1 er LMA allowed region, the comparison of the 
KamLAND and the SNO CC data will determine the total 8 B flux with an uncertainty that is less than 14%; the 
uncertainty can be less than 10% over a significant fraction of the current la allowed domain. About ~ 6% of the 
current estimated uncertainty is due to the experimental error in the SNO CC measurement, which hopefully will be 
reduced as more CC data are accumulated. 

For some purposes, it is convenient to know what is the average expected accuracy in the determination of /b- We 
have computed this average with the aid of a Monte Carlo sampling of the currently allowed solar neutrino oscillation 
parameters shown in Fig. fl\. The average positive and negative uncertainties are approximately equal (cf. Tabic pi) . 
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FIG. 4. Accuracy of determining the total 8 B neutrino flux. The figure displays la contours for the percentage accuracy 
in determining the 8 B flux that can be obtained from the combined SNO CC and KamLAND data. The uncertainties were 
calculated from Eq. ( |l8| ) and the currently allowed regions for the neutrino oscillation parameters were obtained by a global fit 
to all of the allowed solar and reactor data (cf. Fig. fiT). 



We find an average la uncertainty of 

*(/b) = 9.6%. (19) 

A very significant component of ct(/b) comes from cr((P ee )sNO, KamLAND)- We find, averaged over the current la solar 
neutrino oscillation region, 

a((P ee )sNO, KamLAND ) = 7%. (20) 

The sterile neutrino contribution to the 8 B neutrino flux can be determined by subtracting the active neutrino flux 
from the total neutrino flux. Thus, 

, RsSo XP ~ R^ cxp (Pee(Am^ tan 2 0)) SNO 

/B.sterile - i^M^^^^ > W 

Using, as described above, the la errors for the total and the active fluxes of 9.7% and 8% (assumed for the SNO NC 
measurement), respectively, we estimate that the sterile component of the 8 B neutrino flux can be determined to a 
precision of about 12.5%. 



IV. HOW WELL CAN WE DETERMINE THE TOTAL 7 Be FLUX USING CC (RADIOCHEMICAL) 

EXPERIMENTS? 

We show in this section how the total 7 Be solar neutrino flux can be determined, with the judicial aid of other 
neutrino fluxes predicted by the standard solar model [llj], by combining the results of the GALLEX jl5|],GNO |l6| , 
and SAGE Jl3] gallium solar neutrino experiments with the KamLAND and SNO CC measurements. We shall also 
explore the extent to which the Chlorine experiment Jl2 ,Q can provide independent information about the 7 Be solar 
neutrino flux. 

We limit ourselves in this section to detectors that only observe v e or i> e , specifically, we consider here only the 
radiochemical gallium and chlorine experiments and the reactor anti-neutrino experiment, KamLAND. This limitation 
simplifies the calculations with respect to the role of the sterile neutrinos. However, the radiochemical experiments 
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TABLE III. Gallium neutrino capture rates and solar neutrino fluxes. The table presents the predicted standard 
solar model jll| neutrino fluxes and the calculated gallium neutrino capture rates, with la uncertainties from all sources 
(combined quadratically) . The neutrino fluxes are the same as in the original BPOO model pi| . The last column of the table 
presents the capture rate for gallium predicted by the best fit LMA solution. The total rates were calculated using the neutrino 
absorption cross sections and their uncertainties that are given in Ref. jjj. 
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suffer from the disadvantage of a lack of energy discrimination, which introduces uncertainties involving the roles of 
the pp, pep, and CNO neutrinos. 



We begin by describing in Sec. IV A the general procedure for determining the 7 Be neutrino flux. We then evaluate 
in Sec. IV E| the principal sources of error, taking account of experimental and theoretical uncertainties as well a s the 
possibility of an appreciable sterile neutrino component in the incident solar neutrino flux. We present in Sec. IV C 



the numerical results for the uncertainties due to different factors and evaluate the overall accuracy with which the 
total 7 Be flux can be determined. 

Using data from either the gallium or the Chlorine Jl2| experiments, the same proc edure can be applied for inferring 
the 7 Be neutrino flux. For simplicity, we describe the proce dure i n Sec. [TV A[ -Sec. IV C with reference to the more 
promising case provided by the gallium experiments. In Sec. IV D, we investigate how accurately one can determine 
the Be fl ux u sing data from the Chlorine experiment instead of the gallium experiment. 

In Sec. VB, we use the techniques developed in this section to explore the accuracy with which KamLAND and 
BOREXINO can determine f Bc - 



A. Procedure for determining the total 'Be solar neutrino flux 



Table III 



shows the neutrino fluxes and the event rates in the gallium solar neutrino experiments that are predicted 
by the standard solar model |ll],^2| ■ The table also shows the event rate predicted by the best-fit LMA solution. 
From Table III it is clear that one must make a strong assumption about the pp neutrino flux in order to determine 
the 7 Be flux. One must also make assumptions regarding the best-value and the uncertainties in the CNO and pep 
fluxes. This situation is different than the purely empirical procedure described in Sec. HI for determining the 8 B 
neutrino flux; the 8 B solar neutrino flux can be determined independent of all considerations regarding the standard 
solar model. 

We assume throughout this section the correctness of the calculated standard solar model |Tl| values for the 
neutrino fluxes and their uncertainties, except for the Be and 8 B fluxes which we want to determine from solar 
neutrino experiments^]. 

Fortunately, as we shall see, if we assume the SSM predictions and their uncertainties for all but the 7 Be and 8 B 
fluxes, then we can infer an interesting range for the total solar 7 Be neutrino flux if we use data from the gallium 



5 The ultimate astronomical goal of solar neutrino experiments is to determine all of the solar neutrino fluxes directly from 
experiment, but there are too few experimental constraints to make this possible at the present time. 
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tan 2 

FIG. 5. Isocontours for the total 7 Be flux. The figure compares isocontours for the total 7 Be flux assuming oscillations 
between v & and purely active neutrinos (thicker lines) or oscillations between u e and a 75% active-25% sterile admixture 
(thinner lines). The results were obtained by solving Eq. (^3|) as described in the text. The differences between the pure active 
contours and the 75%-25% admixtures are less than 2% within the currently allowed 3<r solution space, which is shown by the 
dotted line in the figure (cf. Fig. [l]). The dotted curve represents the 3cr allowed regions from the analysis of the solar data. 



experiments. The situation is less promising if we use data from the Chlorine experiment rather than the gallium 
experiments. 

Suppose that KamLAND observes a signal that corresponds to V e oscillations with parameters 
(Am 2 , tanf? 2 )i<amLAND, then the expected event rate in the gallium experiments is a sum of the contributions from 
the different neutrino fluxes, namely, 



D £ D s B,SSM /n /a -2 , 2 o\ \ S B , r D 7 Bo,SSM /n) i \ -1 , 2 n\ \ 7 Bc 

-RGa = jB-R Ga (P ee (Am ,tan ) KamLAND / Ga + /Be P Ga (P ee (Am , tan V J KamLAND ) Ga 

+ / ! ^Ga SSM <^ e (A77 l 2 , tan 2 $) KamLAND) Ga' 



(22) 



In the last term in Eq. (|22j), we include the contributions from hep, pep, CNO and pp neutrinos. By analogy with 
Eq. (Q) , we have defined the factors fi as the ratios between the "true" solar neutrino fluxes and the fluxes predicted 
by the standard solar model. 

We can solve Eq. (22) for the 7 Be solar neutrino flux as follows. We substitute into Eq. (|22|) the value of /b 
determined, independent of the solar model, from the KamLAND and SNO CC measurements, as discussed in the 
previous section (Sec. III). We also assume as a first approximation that all the solar neutrino fluxes but the 8 B and 
7 Be fluxes are equal to the values predicted by the SSM; we investigate later the accuracy of this approximation. 
With these assumptions, we can then solve for f-Q e by equating i?Ga = ^Ga = ^.4 ± 4.7 SNU. Thus 



/bc — 
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R, 



7 Bc,SSM 
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P ee (Am 2 , tan 2 6) KamLAND ) G 



7 Bc 
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Z^Ga 



Pee (Am 2 , tan 2 ^KamLAND) G a 



B,SSM ^SNQ XP (-Pee(Am 2 ,tan 2 fl) K amLAND)G" 



R, 



Ga 



Pj?NO SM (Pee (Am 2 , tan 2 



I KamLAND /SNO 



(23) 



Figure || shows that for the 7 Be neutrino flux the uncertainty due to the potential effect of sterile neutrinos is small. 
The figure shows the isocontours of /bc in the LMA region for the case of purely active case (thicker lines) and for 
the mixture of active and sterile neutrinos with cos 2 r\ = 0.75 (thinner lines). Within the la (3<r) parameter region, 
the difference in the value of /b between the two oscillation scenarios is less than 1% (2%). 
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B. Principal sources of uncertainty in determining the 'Be total flux 



The uncertainty in the inferred total 7 Be solar neutrino flux can be estimated from Eq. (p3|). Including just the 
largest contributions, we can write the fractional uncertainty in the total 7 Be neutrino flux as 



c(/Be) 



/Be 



^ — ^(/Bc) G a : cxp + <5(/bc) Ga.crosssection + ^(/Bc) K amLAND + ^(/bc) C NO + 

+ smaller terms (24) 



where <K/Be)G aoxp is the uncertainty from the experimental error of the gallium rate, #(/Be)Ga cross section * s * ne 



theoretical uncertainty from the calculated v e -Ga. absorption cross sections |44f] , S (/Be)KamLAND * s ^ ne uncertainty 
arising from the allowed range of neutrino parameters determined by KamLAND and SNO [which affects all of the 
averaged survival probabilities in Eq. and 8(fB e )cuQ ^ s ^ ne uncertainty due to the quoted errors in the standard 



solar model calculation of the CNO fluxes (see Table III and Ref. [fLl|]). We have omitted from Eq. ([24]) a number 



of sources of error that contribute only relatively small uncertainties. The omitted sources of error [and their range 
of contributed uncertainties for a representative set of oscillation parameters] include the experimental error in the 
SNO CC measurements [(1.2 ± 0.6)%], the uncertainty in the theoretical absorption cross section for the SNO CC 
experiment [(0.7±0.2)%], and the uncertainty in standard solar model calculation of the pp neutrino flux [(0.4±0.1)%]. 

We now discuss how we estimate the uncertainties in /Be- Because they require special treatment, we first discuss 
the uncertainty arising from the theoretical neutrino capture cross sections for gallium and then discuss the uncertainty 
resulting from the range of allowed neutrino parameters determined by KamLAND^. 

To evaluate the uncertainties arising from the gallium absorption cross sections for neutrino sources with continuous 
energy distributions, we use Tables 2-4 of Ref. E|; these tables give the best-estimate and the ±3cr limits for the 
theoretical cross sections. For the neutrino lines from 7 Be and pep, we have checked that the shape of the line E|] 
does not affect significantly the error estimate. Therefore, we use for the 7 Be and pep lines the error estimates given in 
Eq. (41) and Eq. (42) of Ref . Q. We have adopted the conservative procedure described in Sec. XII. A. 4 of Ref. Q, 
in which the uncertainties in all of the low energy (< 2 MeV) cross sections are fully correlated, while the uncertainties 
for the ( 8 B) neutrinos above 2 MeV are treated separately. All of the cross sections for low energy neutrinos move 
up or down together, reflecting the fact that the dominant uncertainties for low energy neutrinos are common to all 
sources. For higher-energy neutrinos, a number of excited states dominate the calculated absorption cross section 
(For a more explicit description of how the cross section errors are treated, see the Appendix.). 

To calculate the uncertainty associated with the range of allowed neutrino parameters determined by KamLAND, 
we first solve Eq. ( p3| ) with the appropriate average survival probabilities computed for each neutrino oscillation 
point (Am 2 , tan 2 9) in the allowed region (cf. Fig. Q). We consider active-sterile neutrino admixtures permitted 
by the currently allowed global oscillation solution [see Eq. fll3|)]. The solution of Eq. (^3|) determines /B e (best fit) 
for that particular point in oscillation parameter space. Then we construct a la allowed region, a set of p oints 



(Am 2 , tan 2 9), a round the chosen point using the simulated characteristics of the KamLAND experiment (see Sec. Ill C 
and Refs. M , gg , B9[ ) . We define 5 ( /bc ) KamLAND for the chosen (Am 2 , tan 2 9) to be the maximum (or minimum) value 
of (/Be(An?, tar?l?) — /Be(best fit)) //B e (best fit). In practice, the inclusion of sterile neutrinos only slightly affects 
the computed range of £>(/Bc)KamLAND- 

For all other quantities, we estimate the associated uncertainty at a particular point (Am 2 , tan 2 9) in the following 
way. Given a source of uncertainty i (for example, the measured capture rate for the gallium experiments) with la 
error Cj, we obtain #(/s e )i from the relation 

<K/Be)i = 7 7~<\ ■ (25) 

/Be l ) 



Here we denote by /bc(*) the value of /Be obtained from Eq. (|23|) when all the parameters are assigned their best- 
estimate values. The quantity /b c (« ± &i) is calculated from Eq. ( p3| ) using the best-estimate values of all variables 
except i; the variable i is set equal to its best-fit value ± the corresponding la uncertainty. In calculating <5(/Be)cNO> 
we shift the three CNO neutrino fluxes by ±la simultaneously and in the same direction, reflecting the correlation 
between the CNO fluxes in the standard solar model. 



6 The SNO CC results are used to select the first quadrant for 8 (cf. Fig. ^j), but for brevity we refer to the range of neutrino 
parameters determined by KamLAND. 
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The uncertainty, <5(/Bc)i, that is calculated from Eq. ( p5| ) will in general depend upon the assumed value 
for (Am 2 , tan 2 9) within the KamLAND allowed region. This dependence persists even if ai is independent of 
(Am 2 , tan 2 9) (which is true, e.g., for the measured gallium capture rate). In fact, the positive and negative val- 
ues for <5(/Be)i will generally not be equal. 

TABLE IV. Best-fit values of fu e and associated uncertainties (gallium based simulation). The table shows the best-fit values 
and associated uncertainties that are obtained by solving Eq. (E3I) for a representative set of oscillation parameters within the 
expected KamLAND and SNO allowed region (cf. Fig. hh. 











Uncertainties <5(/Be)i (%) 
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C. The accuracy with which the total 7 Be flux can be determined 



Table IV presents the calculated uncertainties and the best-fit values of /b g for a representative set of possible 
neutrino oscillation parameters, (Am 2 , tan 2 9), that may be obtained from the KamLAND measurements . The 
largest uncertainty, ~ 22%, is due to the experimental error on the measured gallium rate p4]-^6|. The two next 
largest uncertainties, both ~ 12%, arise from the theoretical calculation of the gallium absorption cross sections and 
the simulated errors in the KamLAND measurements. The rather large uncertainty due to the gallium cross sections 
requires explanation since the uncertainties on the individual cross sections are much smaller p4| [e.g., 2.3% for pp 
neutrinos and 5% (average) for 7 Be neutrinos]. The amplification in the error due to the cross sections arises because 
all of the low energy cross section errors are fully correlated. The gallium cross section errors are added linearly in 



calculating SfBe [cf. Eq. ( |23j )] rather than being combined quadratically. 

We have also computed a representative (average) error in the determination of the total 7 Be flux by a Monte Carlo 
sampling of the allowed KamLAND oscillation region shown in Fig. [j]. We find 

/bc - /sew fit [1.00ig;g] ■ (26) 

Figure || shows the contours of the maximum percentage deviation (in absolute value) from the local best-fit value of 
/bc- The figure shows that within the currently allowed la solar neutrino oscillation region the expected uncertainty 
in the determination of the 7 Be flux is of the order of 25% — 35% [in agreement with Eq. (pF 



D. Can one use the Chlorine experiment to determine the 7 Be solar neutrino flux? 



We can derive an expression for /b g in terms of the measured event rate in the Chlorine |12| , [l3f solar experiment. 
Replacing "gallium" by "chlorine" everywhere in Sec. IV A -Sec. IV C, we find 
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FIG. 6. Percentage error in determining the total 7 Be flux. The figure shows two conto urs (in % ) for the uncertainty in 
determining the total 7 Be flux. The uncertainties were calculated as described in Sec. |[V A and Sec. IV B, using a combined 
analysis of SNO CC and gallium data together with simulated KamLAND data. The curves labeled la, 2a, and 3a represent 
allowed regions from a global analysis of the available solar and reactor data (see Fig. |l|). 
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(27) 



where i?^ p = 



2.56 ± 0.23 SNU |12|| . Equation (|27| ) is the analogue for the Chlorine experiment of the previously 
derived Eq. ( p3| ) that was used in the discussion of extracting /B e for the gallium experiments. 

Table |V lists the best-fit values of /b c that were obtained by solving Eq. (^) for different pairs of oscillation 



parameters, [Am , tan 



The best-fit values of /b 6 obtained from the the Chlorine experiment are smaller than 



the best-fit values inferred using the gallium data (cf . Table |y| and Table 0) . Using the chlorine data, one can even 
get negative (i.e., unphysical) solutions for /b c . 

The basic reason for the difficulty in determining /b c is that the large expected rate from the 8 B neutrino flux, 
even after reduction due to neutrino oscillations, can account for all the observed rate in the Chlorine experiment p6[ . 
The contribution from the 7 Be neutrinos is obtained by subtracting the large and somewhat uncertain expected 
contribution from 8 B neutrinos (and the contributions from CNO and pep neutrinos that are expected to be much 
less important) from the total measure chlorine rate. The result is a rather small and uncertain remainder, which can 
be attributed to 7 Be neutrinos. 

Table ^ also presents the most important sources of uncertainty for inferring the value of /bc using the chlorine, 
rather than the gallium, data. Since the inferred best-fit values for /b c that are obtained using the chlorine data are 
very small (or even negative), we show in Table [V] the associated shifts in the prediction of /b c - For the much more 
reliable inferences from gallium data, we show instead in Table |^ the percentage shifts, not the actual numerical 
shifts. 

The largest uncertainties in determining /b c using the chlorine data are caused by the experimental errors in the 
chlorine absorption rate (~ 0.35) and the SNO CC absorption rate (~ 0.25). The uncertainty in the chlorine neutrino 
absorption cross sections are also significant (0.03 to 0.35, depending upon the neutrino oscillation parameters). To 
be explicit, the uncertainties in /b c due to the theoretical uncertainties in calculating the chlorine neutrino absorption 
cross sections that were used in constructing Table [v] were evaluated from the following equation, 



A(/Be)ci,c. S . 
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=pep, 7 Be,CNO 
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(28) 



16 



TABLE V. Best-fit values of /b c and associated uncertainties (chlorine based simulation). This table is similar to Table [lV| 
except that for Table ^ data from the Chlorine experiment were used instead of data from the gallium experiment. Because 
the inferred values of /b c are small and very uncertain using data from the Chlorine experiment, Table VI p resents uncertainties 
as the numerical shift in the best-fit values (not as percentage uncertainties, which are given in Table Evj). 



A - 2 

Am 


2 n 

tan 8 


/Be 






Uncertainties A(/e e )i 






CI, exp 


CI, c.s. 


SNO,exp 


SNO,c.s. 


CNO, flux 


Total 


5.0 x 10" 5 


4.2 x 10 _1 


0.17 


0.37 


0.25 


0.26 


0.07 


0.08 


0.50 


5.0 x 10" 5 


5.01 x 10" 1 


0.20 


0.38 


0.29 


0.21 


0.07 


0.08 


0.53 


5.0 x 10" 5 


2.51 x 10" 1 


0.07 


0.32 


0.10 


0.18 


0.06 


0.08 


0.39 


7.94 x 10~ 5 


4.2 x 10 _1 


0.12 


0.36 


0.17 


0.20 


0.07 


0.08 


0.46 


7.94 x 10~ 5 


5.01 x 10 _1 


0.16 


0.38 


0.22 


0.21 


0.07 


0.08 


0.49 


7.94 x 10~ 5 


2.51 x 10 _1 


-0.10 


0.31 


0.10 


0.17 


0.06 


0.08 


0.37 


3.16 x 10~ 5 


4.2 x 10 _1 


0.25 


0.39 


0.34 


0.21 


0.07 


0.07 


0.57 


3.16 x 10~ 5 


5.01 x 10" 1 


0.26 


0.40 


0.37 


0.22 


0.07 


0.07 


0.60 


3.16 x 10 -5 


2.51 x 10" 1 


0.16 


0.36 


0.25 


0.18 


0.06 


0.07 


0.46 



(See the Appendix for more details regarding the treatment of the cross section errors.) Even the CNO fluxes 
(~ 0.08) and the SNO CC absorption cross section (~ 0.10) contribute non-negligible errors. The uncertainty from 
sterile neutrinos is small (~ 0.04) and does not significantly affect the total uncertainty. 

Remarkably, the individual uncertainties in /b c which arise from several different sources are larger than the current 
best- fit values of /bc (see Table [v]) . 

For the chlorine based determination, we have computed a representative (average) shift in f^ e by a Monte Carlo 
sampling of the allowed KamLAND oscillation region shown in Fig. [|. We find 

A/bc = 0.49, la. (29) 

This uncertainty is so large as to render not very useful the determination of /b c with the aid of chlorine data. 



V. HOW WELL CAN KamLAND PLUS v - e SCATTERING EXPERIMENTS DETERMINE THE 

TOTAL 8 B AND 7 Be FLUXES? 



In this section, we show how data from KamLAND can be combined with v — e scattering data obtained with 
the Super-Ka mioka nde and BOREXINO experiments in order to determine, respectively, the total 8 B (Sec. V A ) 
and 7 Be (Sec. VB) solar neutrino fluxes. As discussed before, the advantage of using purely CC measurements to 
determine the fluxes is that the answers depend only mildly on the active-sterile admixture. On the other hand, 
v — e scattering measurements have the advantage of smaller uncertainties in the interaction cross sections. Moreover, 
for the determination of the 7 Be flux, BOREXINO depends less strongly on other neutrino fluxes predicted by the 
standard solar model than do the radiochemical experiments, Chlorine, GALLEX, SAGE, and GNO. 



A. How well can KamLAND and Super— Kamiokande determine the total 8 B flux? 

In this section, we show how data from the Super-Kamiokande and KamLAND experiments can be combined to 
measure the total 8 B neutrino flux. We shall see that the determination using Super-Kamiokande and KamLAND 
yields, on average, a value for the total flux than is comparable in precision to what is expected to be obtained using 
SNO and KamLAND. The systematic uncertainties are different in the two experiments, SNO and Super-Kamiokande. 
Therefore, it will be important to compare the total 8 B neutrino flux that is inferred using Super-Kamiokande and 
KamLAND with the value that is obtained using SNO and KamLAND. If the uncertainties are correctly estimated, 
then the two methods should agree within their quoted errors. 

The advantage of using purely charged current measurements with SNO and KamLAND to determine the total 
8 B flux is that the answer depends only very slightly upon the unknown active-sterile mixture (as discussed in 
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Sec. [II). The principal advantages of the Super- Kamiokande experiment in the present context is that the neutrino 
interaction cross section is accurately known |47| and the statistical and systematic errors have already been extensively 
investigated ||. However, as we shall see below, when using Super-Kamiokande there is a significant uncertainty (+7% 
and —0%) in the total 8 B flux due to the active-sterile mixture, at least with our present knowledge of r\. 
The average survival probability for 8 B solar neutrinos can be written in the form 

(P(Am 2 , tan 2 9, 77)KamLAND)supcr-Kamiokandc = (Pee ( Am 2 , tan 2 9, ^KamLAND) Super— Kamiokande 

+ (?"-Pex(Am 2 , tan 2 0, 77)KamLAND)supcr-Kamiokandc, (30) 

where P is the probability of oscillating into active neutrinos and r = a^/a e ~ 0.15 is the ratio of the the v e — e 
and Ufx — 6 elastic scattering cross-sections. The expected sensitivity of /b to the principal sources of errors can be 
calculated from the following equation: 



c(/b) \ | °"(^Supcr-Kamiokandc) \ / &((P) Super-Kamiokande, KamLAND 



/B J 

\ "^Supcr— Kamiokande / V (^ > ) Super — Kamiokande, KamLAND ) 



(31) 



We suppose that KamLAND will observe (with associated uncertainties that we simulate with a Monte Carlo) the 
rate predicted for the global best-fit point shown in Fig. [I]. For purely active oscillations, we find that 

/aetive B - 1-07 [lt° f 2 l = 1 [1 ± 0.054] , COS 2 T, = 1.0, (32) 

where the first error corresponds to the Super-Kamiokande experimental uncertainty and the second error is caused 
by the finite size of the allowed KamLAND region in oscillation parameter space. The combined KamLAND and 
Super-Kamiokande measurement will, for purely active neutrinos, yield a determination for the total 8 B flux that is 
more accurate than can be obtained with the SNO CC measurement and KamLAND. Within the 1 a LMA region, 
the average uncertainty in the combined KamLAND and Super-Kamiokande measurement for the total active 8 B 
neutrino flux is expected to be 

^Cfaetive fi) = 0.06. (33) 

For purely active neutrinos, Super-Kamiokande and KamLAND may provide us with a more accurate determination 
of the total 8 B neutrino flux than SNO CC and KamLAND [cf. Eq. © and Eq. ©]. 

Allowing for the currently allowed la (25%) sterile admixture, we find for the best fit point 

f — 1 07 Pi +0-037 +0.040 +0.07] ^J^xnTr (o A \ 

/total B - 1-07 [1_ .028 -0.046 -O.OoJ ) COS Tj > U.75 (d4) 

The last error in Eq. ( p4] ) corresponds to the present la uncertainty from the active-sterile admixture cos 2 77 > 0.75. 
The average uncertainty in the combined KamLAND and Super-Kamiokande measurement for /total b is expected to 
be ig%, within the 1 a LMA region and if the sterile admixture is as large as currently allowed. 



B. How well can KamLAND and BOREXINO determine the total Be flux? 



In this section, we show how data from the KamLAND and BOREXINO experiments can be combined to measure 
the total 7 Be neutrino flux. The principal advantage of using the BOREXINO experiment for this purpose is that the 
signal in the BOREXINO experiment is predicted to be dominated by 7 Be neutrinos, whereas 7 Be solar neutrinos 
are expected to contribute only a relatively small (and unlabeled) fraction to the observed event rate in the gallium 
and chlorine radiochemical measurements (see Table |n|)Q 

However, unlike the cases involving the gallium and chlorine experiments that were discussed in Sec. [PV| , which 
include only i/ e (CC) absorption, the BOREXINO experiment detects both i/ e — e scattering and — e and v r — e 
scattering. One must consider in the present case the extent to which the active-sterile neutrino admixture influences 
the detected event rate for each set of oscillation parameters, ( Atti 2 , tan w )kland 1 determined by KamLAND. As we 



7 The BOREXINO detector can measure the energy of the recoil electrons produced by v — e scattering. The radiochemical 
detectors do not have energy resolution, only an energy threshold. 
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shall see quantitatively in the following discussion, this uncertainty regarding the sterile admixture does not decrease 
significantly the overall accuracy of the inferred total Be neutrino flux. 

The average survival probability for active neutrinos can be written conveniently in the form 

(P(Am 2 , tan 2 8, 7?) KamLAND ) BOREXINO = (-Pee(Am 2 , tan 2 0, ff) KamL AND ) BOREXINO 

+ (rP ex (Am 2 ,tan 2 0,77)KamLAND) B OREXINO' ■ ( 35 ) 

where P^ >i is the oscillation probability of neutrino fluxes of source 0j into active neutrinos and r = a^j 1 <j e ~ 0.15 is 
the ratio of the the v e — e and — e elastic scattering cross-sections. 

The expression for /s e has the same form for the case in which the KamLAND and BOREXINO experiments are 
considered together as for thcpreviously discussed cases involving the gallium experiments [see Eq. (p3)] and the 
chlorine experiment [see Eq. (P7|) l . Using the average survival probability defined in Eq. ( |35"1) , we can write 



f 1 

/Be 



p 7 Bo,SSM //DA _ 2 + 2 n \ \ 7 Be 
BOREXINO \V ' tan ' ^/KamLAND/ BOREXINO 



noxp 

-"•BOREXINO 



~ ^BOREXINO (-P(Am 2 , tan 2 9, rf) KamLAND ) B OREXINO ( 36 ) 

i 

In the last term of Eq. ([36]), we include the contributions from pp, and CNO neutrinos. The contribution from 8 B 
neutrinos is negligible because the observed 8 B neutrino flux is about a factor of 10 3 smaller than the predicted SSM 
7 Be neutrino flux and because 8 B neutrinos primarily produce high energy recoils electrons that the BOREXINO 
detector can distinguish from the low energy recoil electrons produced by 7 Be neutrinos. 
The expected sensitivity of /Be to different sources of errors is given by 



Q"(/bc) 

/Be 



BOREXINO, exp + S(f Bc ) 

KamLAND + ^(/Be)cNO + smaller terms. (37) 



In order to estimate the accuracy of this method for determining /b c , we suppose that KamLAND will observe 
(with associated uncertainties that we simulate with a Monte Carlo) the rate predicted for the global best-fit point^j 
shown in Fig. [l] and that BOREXINO also will observe a signal corresponding to this best-fit point (with associated 
uncertainties). In the absence of any published data, we estimate that-based upon experience in previous solar 
neutrino experiments P,^| Jl^Jl-| |l6|-that BOREXINO will achieve a systematic uncertainty of between 5% and 10% 
. We estimate in this way that the combined experiments will yield a determination of /b 6 that is, for purely active 
oscillations, 

/active Be = 1-00 [1 ± 0.05(0.1) ± 0.020 ± 0.020] = 1 [1 ± 0.057(±0.103)] , cos 2 j) = 1.0, (38) 

where the first error corresponds to the BOREXINO experimental uncertainty, the second to the uncertainty in the 
reconstructed KamLAND region, and the last error is due to the theoretical uncertainty in the predicted CNO fluxes. 
Within the 1 a LMA region, the average uncertainty in the combined KamLAND and BOREXINO measurement for 
the total active 7 Be neutrino flux is expected to be 

^/active Be) = 0.06(0.105) . (39) 

If we allow for a 25% sterile admixture, we find 

/total Be = LOO [1 ± 0.05(0.1) ± 0.02 ± 0.020 ±°;°j>] . (40) 

The last error in Eq. ( [hil ) corresponds to the present la uncertainty from the active-sterile admixture cos 2 r\ > 0.75. 
The numbers in parentheses in Eq. (|3^) and Eq. ( f40| ) correspond to assuming the larger systematic uncertainty, ±10%, 
for the BOREXINO measured rate. 

For /total Be, the average uncertainty in the combined KamLAND and BOREXINO measurement is expected to be 
tl{tll)%, within the 1 a LMA region and if the sterile admixture is as large as currently allowed, 



8 The predicted rate in units of the expected SSM rate is J?borexino = 0.64, see Ref. |36| 
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VI. DETERMINATION OF THE pp NEUTRINO FLUX 



In this section, we analyze three strategies for determining the total pp solar n eutrin o flux without requiring a 
dedicated experiment that measures only the pp neutrinos. We first describe in Sec. VI A how one can make a crude 
determination of the pp neutrino flux using the measured Gallium, Chlorine, and SNO (CC) event rates f3|, |l2] , [l5|[L6[ 
and the standard solar model predictions for all but the 8 B, 7 Be, and pp solar neutrino fluxes. This part of the 
discussion is similar to the analysis described in Ref. || , although we evaluate explicitly the uncertainty caused by the 
finite size of the allowed region in oscillation parameter space . We then determine in Sec. VI B how well one can infer 
the pp flux using just the Gallium and the SNO measurements and the BPOO predictions for the other neutrino fluxes, 
especially the Be neutrin o flux . The u nknow n sterile-active mixtur e con tributes only a negligible error using the 



strategies described in Sec. VIA and Sec. VI B. Next, we show in Sec. VIC how the precision of the determination of 
the pp flux can be improved by using, in the future, the results from KamLAND (to constrain the neutrino oscillation 
parameters) and from BOREXINO (to constrain the 7 Be neutrino flux). The principle of this strategy has also been 
discussed by the SAGE collaboration || . 

The theoretical error on the pp neutrino flux is ±1% (see Ref. [0), which is an order of magnitude smaller than 
the estimated errors that we find in this section on the empirically derived values of f pp . To achieve a precision of 
10% or better in the determination of the solar neutrino pp flux will require a dedicated and accurate experiment 
devoted to measuring the pp flux. 

Throughout this section, we treat the pp flux and the pep flux as a single variable because they are so closely linked 
physically [ fl8f . The two reactions are linked because the pep reaction is obtained from the pp reaction by exchanging 
a positron in the final state with an electron in the initial state. Thus the rates for the two reactions are proportional 
to each other to high accuracy p9[. For convenience, we denote the sum of the pp and pep fluxes as simply pp. 



A. Using Chlorine, Gallium, and SNO data and BPOO predictions 



In this subsection, we show how data from the chlorine, gallium, and SNO experiments can be combined with the 
BPOO predictions for the CNO and hep fluxes to determine the total pp neutrino flux. 

The reduced pp solar neutrino flux, f pp defined (by analogy with /b) with respect to the predicted standard solar 
model sum of the pp and pep fluxes, can be written as: 



fp 
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^ a SSM (P ee (Am2,tan^))^ a 
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(41) 



Here the sum over i in Eq. (41) refers to the three CNO neutrino sources and the hep neutrinos (see Table HI). 

We insert in Eq. (41) the expression for /bc in terms of the Chlorine and SNO rates and the standard solar model 
CNO and hep neutrino fluxes. Explicitly, 



/Be 
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3c,SSM 
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For the special case of the best fit point in the LMA solution region, we find 

f - 1 41fl± 08 +0 03+0 009 ± 06 +006 +0 04s ) 

Jpp — A-H-m =c U.UO _o.05 -0.007 ^ u - uu -0.19 -0.il> 



(42) 



(43) 



The first error in Eq. ([43]), (±0.08), results from the weighted average experimental error for the Gallium exper- 
iment HHHl. The second error, (^q'q 3 ), reflects the uncertainties in the calculated neutrino absorption cross 
sections on gallium jH. The third error, (ig'oo?)' ^ s causecl by the uncertainties in the calculated standard solar 
model CNO fluxes |11[|. The fourth error, (±0.06), derives from the measurement errors in the SNO CC experimental 
data Q and neutrino absorption cross section [fHf . The fifth error, (j^ ^g), results from the experimental error for the 



Chlorine event rate | 
cross sections for Ch 



12| and the sixth error, (io?i)i reflects the uncertainties in the calculated neutrino absorption 
orine The physical constraint, /b c > 0, was imposed in obtaining these errors. 
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Within the la LMA allowed region, we find that the central value of f pp can vary in the range 1.32 < f pp < 1.50 
and the best-fit value can be determined with an average uncertainty of j^4%- The range of central values of f pp 
always exceeds the BP00 value because the observed Chlorine rate implies a low value of /v Bc . Since the largest 
contributions to the Gallium rate are from pp neutrino and 7 Be neutrinos, a low 7 Be flux implies a relatively high pp 
flux. 

We can summarize as follows the results of the simulations carried out within the la allowed LMA region for the 
pp flux: 

f pp = 1.41(1 ±0.06^24) = 1-41(1^). (44) 

The first error in Eq. ( f44| ) results from the allowed range of neutrino oscillation parameters and the second error 
results from the uncertainties in all other recognized sources of errors, combined quadratically. The result given in 
Eq. ([wj) should be compared with the quoted estimate by the SAGE collaboration ||, f pp = 1.29(1 ± 0.23). This 
agreement is very welcome since the SAGE paper points out that they made "Several approximations... whose nature 
cannot be easily quantified." 



B. Using Gallium and SNO data and BPOO predictions 

In this subsection, we determine the range of allowed values for f pp using the BPOO predictions (and uncertainties) 
for the 7 Be neutrinos as well as the CNO and hep neutrinos. We use in this subsection data only from the Gallium 
and SNO solar neutrino experiments. 



Following the same line of reasoning as in section VI A , wc find for the best fit point in the LMA allowed oscillation 
region, 

fpp = 1-05(1 ± 0.11 ±g;g| t Q a °ol ± 0.0071°;°*). (45) 



Just as in Eq. (E3), the first error in Eq. (f45|), (±0.11), is the experimental error from the weighted average event rate 

K).05\ 
-0.08J 



in the gallium experiments, the second error, (iglos)) * s ^ ue t° ^ ne neutrino absorption cross section on gallium, the 



third error, (i°!o3 )> ^ s due *° the uncertainties in the BPOO predictions of the CNO neutrino fluxes, the fourth error, 
(±0.007), contains the uncertainty in the SNO CC experimental rate and the calculated absorption cross sections on 
deuterium. The uncertainty in the BPOO prediction for the 7 Be neutrino flux contributes the last error, (i° 02)1 m 
Eq. @. 

Within the la allowed LMA region, we find the central value of f pp varies in the range 0.93 < f pp < 1.16; the 
best- fit value of f pp can be determined with an average uncertainty of ±14%. We can summarize the determination 
of f pp as follows: 

/pp = 1.05(lt°;if ± °- 14 ) = l-05(! ± °- 18 )> ( 46 ) 

where the first error is the uncertainty due to the neutrino oscillation parameters and the second error contains the 
uncertainties from all other sources. 



C. Using BOREXINO, KamLAND, Gallium and SNO data and BPOO predictions 

In this subsection, we show how the uncertainty in / pp could be improved by using the KamLAND data to determine 
the neutrino oscillation parameters and the BOREXINO data, together with the KamLAND oscillation parameters, 
to constrain /7 Bc . We start with Eq. ( fll| ) but now use /b c as determined from BOREXINO and KamLAND data 
[Eq. ©]. 

We find that for the best-fit LMA solution (assuming that BOREXINO finds the rate expected for this best fit 
point), 

fpp = 1.05(1 ± 0.11 t°;°J ± 0.007 ± 0.05 ± 0.02(0.04) ±°;°°) (47) 

The first error shown in Eq. (p7|), (±0.11), is the experimental error on the measured Gallium event rate, the second 
error, (Iq ,^), represents the uncertainties from the calculated gallium cross sections. The third error, (1q 02)' ' s 
from the predicted CNO fluxes and the fourth error, (±0.007), contains the uncertainties from the SNO experimental 
data and calculated deuterium absorption cross sections. The range of oscillation parameters within the KamLAND 
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reconstructed region gives the fifth error, (±0.05), and the sixth error, (±0.02 [±0.04]), results from the uncertainty in 
the measured BOREXINO event rate [which we take to be 5% (10%)]. The unknown active-sterile admixture (which 
goes in the direction of lowering f pp ) contributes the last error shown in Eq. ([47|). 

On average, the the precision expected using the Gallium, KamLAND, and SNO experiments is 

fpp = 1.05[1.0 + 0.14(0.15)]. (48) 

The dominant sources of error in Eq. (|4S| ) are the experimental error in the gallium rate and the uncertainties in the 
calculated gallium absorption cross section. 



VII. DISCUSSION AND CONCLUSIONS 



In this section, we review and discuss our principal conclusions. We begin in Sec. VII A by summarizing the 
experimentally-allowed range that currently exists for the total 8 B solar neutrino flux, taking account of the possibil- 



ity that sterile neutrinos exist. We then summarize in Sec. VII B how well the total B neutrino flux, and separately 



the sterile B neutrino flux, can be obtained by combining KamLAND and SNO measurements. Next we describe 



in Sec. VII C how well the total 7 Be flux can be determined using data from the KamLAND, gallium (GALLEX, 



SAGE, and GNO), and SNO experiments, i.e., using only CC disappearance experiments. In this same subsection, we 
summarize the measurement accuracy that can be obtained using data from just the KamLAND and the BOREX- 
INO (y — e scattering) experiments. We describe three strategies (making use of existing and future solar neutrino 
measurements and guidance from the BP00 solar model) for determining the pp s olar ne utrino flux without the help 
of an experiment that measures separately the pp neutrinos. We point out in Sec. VII E (and in the Appendix) that 
in order to determine the total 8 B or 7 Be neutrino flux, and a fortiori to determine the sterile contributions to these 
fluxes, the correlations among theoretical errors for neutrino absorpt ion cr oss sections and for solar neutrino fluxes 
must be treated more accurately than in previous discussions. Section VII F summarizes the main focus of the present 
paper. 

We concentrate on procedures to determine experimentally the total solar 8 B neutrino flux and the total solar 7 Be 
neutrino flux in a universe in which sterile neutrinos might exist. Our methods can work only if the LMA solution of 
the solar neutrino problem is correct. 

The numerical values given here for the precision with which different quantities can be determined are obtained 
using simulations of how well different experiments may perform. Therefore, the numerical values are intended only 
as illustrative guides to what may be possible. 



A. The currently allowed B solar neutrino flux if sterile neutrinos exist 



The combined SNO CC data and the Super-Kamiokande v e — e scattering data together yield a widely acclaimed 
agreement between the 8 B solar neutrino flux that is predicted by the standard solar model and the empirically 
inferred 8 B neutrino flux. However, this agreement between solar model prediction and solar neutrino measurement 
is not a unique interpretation of the existing measurements if one allows for the possibility that the incident solar 
neutrino flux could contain a significant component of sterile neutrinos. We show in Table | and in Sec. H that if one 
takes account of the possibility that sterile neutrinos may exist then the total solar 8 B neutrino flux could be as large 
as 2.3 times the flux predicted by the standard solar model. In principle, the existing solar neutrino data could be 
inconsistent with the standard solar model predictions. 



B. Determining the total B solar neutrino flux including sterile neutrinos 

The total 8 B solar neutrino flux, active plus sterile neutrinos, can be determined with a typical la accuracy 
of about 10% by comparin g th e charged current measurements from the KamLAND reactor experiment and the 



SNO experiments (see Sec. III). The active B neutrino flux has been measured by comparing the SNO CC flux 
(which measures v e ) with the Super-Kamiokande neutrino-electron scattering rate (which measures u e plus, with less 
sensitivity, + v T ). The SNO neutral current measurement will provide an additional and, ultimately, more accurate 
measurement of the total active 8 B solar neutrino flux. 

By subtracting the independently-measured active 8 B neutrino flux from the total (active plus sterile) 8 B neutrino 
flux, one can determine empirically the sterile component of the solar neutrino flux. The active 8 B neutrino flux will 
eventually be determined accurately by the SNO neutral current measurement [ELriJ . We estimate therefore that the 
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procedure described here has the potential of measuring, or setting an upper limit on, the sterile component of the 
8 B neutrino flux that is as small as 12% of the total 8 B solar neutrino flux. 

The measurement of the total 8 B neutrino flux, and the sterile component of this flux, are independent of solar 
model considerations. In order to establish the quantitative conclusions, we have performed detailed simulations of the 
accuracy of the KamLAND reactor experiment in determining neutrino oscillation parameters (see Fig. ^) and have 
evaluated the theoretical and experimental uncertainties that affect the different flux determinations (see Sec. Ill C 
and the Appendix). 

The combined measurements of the Super-Kamiokande and KamLAND experiments can be used to determine 
independently a value for the total 8 B neutrino flux. This determination may be as accurate as 6% for purely active 
neutrinos. With the current limits on the active-sterile admixture, the total 8 B neutrino flux could be inferred to 



an accuracy of 9% or better, as described in Sec. VA. It will be important to compare the value of the total 8 B 



neutrino flux inferred by combining the KamLAND and SNO charged current measurements with the value obtained 
using the KamLAND and Super-Kamiokande experiments. This comparison will be an important test of whether the 
systematic uncertainties in the experiment and in the analyses arc understood. 



C. Determining the total 7 Be solar neutrino flux including sterile neutrinos 

The total 7 Be solar neutrino flux, active plus sterile neutrinos, can be determined to a la accuracy of about 30% 



by combining measurements from KamLAND, SNO, and the gallium experiments (see Sec. IV). Unlike the purely 
empirical determination that is possible for the 8 B flux, the measurement of the total 7 Be solar neutrino flux requires 
some assumption regarding the CNO solar neutrino fluxes. In our estimates, we have assumed that the standard solar 



model predictions for the CNO fluxes, and their uncertainties, are at least approximately valid. Table IV shows that 
as long as the CNO fluxes are not a factor of three or more larger than the standard solar model predictions, then 
they will not significantly limit the accuracy with which the total 7 Be neutrino flux can be determined. The measured 
capture rate in the gallium experiments (GALLEX, SAGE, and GNO) currently constitutes the largest recognized 
uncertainty in the determination of the total 7 Be flux by the method described here (see Table IV). The constraints 



provided by the Chlorine experiment are not very useful in providing an accurate determination of the 7 Be neutrino 
flux [see Table and Eq. @]. 

One can also determine the allowed range of the total 7 Be solar neutrino flux using the data from the KamLAND 



reactor experiment and the BOREXINO solar neutrino experiment. We show in Sec. VB that with this method one 
may hope to obtain a la accuracy of 11% or better for the total 7 Be solar neutrino flux, which is more accurate than 
we estimate will be possible with the gallium and chlorine radiochemical experiments. 

One can determine an upper limit to the sterile component of the 7 Be solar neutrino flux by combining the measured 
rate in the neutrino-electron scattering experiment, BOREXINO, with the 7 Be total flux inferred from measurements 
of the KamLAND, SNO, and gallium experiments. If one assumes that the entire signal measured in BOREXINO 
is due to v e , then one obtains a minimum value for the active component of the 7 Be neutrino flux. Subtracting this 
minimum value from the total Be flux, one will obtain an upper limit to the sterile component of the flux. It seems 
unlikely that the procedure described here has the sensitivity to measure a value for the sterile component unless the 
sterile flux is larger than 30% of the total 7 Be neutrino flux. However, limits on the sterile neutrino admixture can 



be obtained from the analysis of B and KamLAND neutrino measurements described in Sec. VII B . 

In order to make a direct and precise measurement of the sterile component of the 7 Be solar neutrino flux, we need 
a charged current measurement of the 7 Be flux. A 7 Be solar neutrino absorption experiment, e.g., with a lithium 
target |32| or with LENS [[53), would make possible an accurate determination of the sterile 7 Be neutrino flux by 
providing a set of experimental constraints that is analogous to what will exist for the SNO, Super-Kamiokande, and 
KamLAND experiments. 



D. Determining the pp neutrino flux 

The theoretical uncertainty in the calculated pp solar neutrino flux is estimated to be only 1% [jll|. Therefore, a 
precise determination of the pp solar neutrino flux will be of great interest as a crucial test of the theory of stellar 
evolution. The measurement of the pp neutrino flux will also provide a critical test of whether the neutrino oscillation 
theory, which works well at energies above 5 MeV, also describes accurately the lower energy neutrino phenomena 



(energies less than 0.4 MeV). We show in Sec. VI that the pp flux can be determined with modest accuracy, of order 
18% to 20%, using a combination of existing experimental data and some guidance from the BP00 standard solar 
model. In the future, it should be possible to determine the pp neutrino flux to an accuracy of 15% using experimental 
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data from the BOREXINO, KamLAND, and SNO solar neutrino experiments and predictions (in a non-critical way) 
from the standard solar model (see Eq. fi^). The unknown flux of sterile neutrinos does not significantly affect the 
quoted estimates on the accuracy with which the pp flux can be determined. To measure the pp flux with an accuracy 
sufficient to test stringently the standard solar model prediction will require a dedicated and accurate experiment 
that measures separately the pp neutrino flux. 

E. Correlation of errors, especially for neutrino absorption cross sections 

In order to determine the total 8 B or 7 Be solar neutrino flux, and their sterile components, one must evaluate 
carefully all known sources of error. In the course of this investigation, we realized that the evaluations of the 
neutrino absorption cross section uncertainties in previous neutrino oscillation studies, including our own, have not 
properly taken account of the correlations among the theoretical uncertainties in the cross section calculations (for 
an insightful discussion of this point, see Ref. p4j). We discuss in the Appendix how the cross section errors can be 
treated more correctly. We also emphasize here that it is necessary to treat as fully correlated the uncertainties in 
the principal CNO neutrino fluxes that are obtained from standard solar model predictions. These effects are small 
compared to other uncertainties in determining solar neutrino oscillation parameters (the principal goal of nearly all 
previous oscillation studies). The correlations among cross section uncertainties become important only when one 
wants to make accurate inferences regarding the neutrino fluxes themselves. 

F. Focus: total fluxes and sterile neutrino fluxes 

The main focus in this paper is on determining experimentally the total 8 B, 7 Be, and pp solar neutrino fluxes in 
order to make possible more precise tests of solar model predictions. We have also shown that the contribution of 
sterile neutrinos to the total flux can be measured, or a useful upper limit can be set, for 8 B solar neutrinos and for 
7 Be solar neutrinos. 
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APPENDIX: 

We determine, as is conventional for many analyses of solar neutrino data, the allowed regions in the neutrino 
oscillation space using a % 2 function that includes all the relevant data. In the construction of the % 2 function, we 
have followed closely the prescription of Ref. |55| (see also Ref. |Q), but we have included some modifications to 
this prescription in order to account in more detail for the energy dependence and the correlation of the cross section 
errors for the Chlorine and gallium solar neutrino experiments. 

For a given experiment j (for example, gallium or chlorine), the expected number of events can be written as 

8 8 
i— 1 k i—1 

where i = 1, 8 labels the solar neutrino fluxes, 4>ii an d k labels the energy bins of energy E^. The quantity Cj(Ek) is 
the cross section for the interaction of a neutrino of energy in the experiment j; P ee ,i(i?fe, Am 2 , 9) is the survival 
probability of v e for a given energy Ek at a specified point, Am 2 , 6, in neutrino oscillation spaceQ. 



The survival probability depends upon the neutrino species, i, because different neutrino species have different probability 
distributions for the location of their production within the Sun. For matter oscillations, the survival probability obviously 
depends upon where in the Sun the neutrino was produced. 
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The error matrix for the neutrino absorption cross sections can be derived using Eq. (Al). Since the errors are 
uncorrelated between different experiments, we need to evaluate the following expression for a specified experiment: 

a 2 (c.s.) = ((R-R°[ 



(A2) 

where R is best-estimate for the (gallium or chlorine) experimental capture rate and the brackets indicate an average 
over the probability distribution of uncertainties in the neutrino cross sections. Writing out the various terms in 
Eq. ([A|), we find 



(7 2 (c.s.) = J2 <Pi( E ki)<P j (E k 2)Pi(E kl )P j (E k2 ) < AC(E kl )AC(E k2 ) > 



(A3) 



In the standard Ref. 



i,j k\,k2 

the cross section error matrix for each experiment is given as: 

^(c.s.)^ = 5 lhl2 j2Y. R ^ R ^ AhiC ii AhiC a 

il i2 

= ^^(AlnCf) 2 



(A4) 



where A In C j is the average uncertainty in th e interaction cross section of neutrinos of flux type i in experiment j. 

Equation (A4) is obtained from Eq. (Al) by neglecting the correlations of the cross section errors among the 
different neutrino fluxes and by neglecting the energy dependence of the cross section errors. For the gallium detector, 
neither of these assumptions is correct and even for chlorine the cross section errors for different neutrino sources are 
strongly correlated. 

For all but the 8 B and hep neutrino fluxes, either ground state to ground state transitions are the only energetically 
possible transitions (which is the case for the chlorine detector) or the ground state to ground state transitions 
dominate (which is the case for the gallium detector) . Thus all the cross sections for the lower energy neutrinos move 
up or down together, proportional to the square of the dominant matrix element. 

When considering the gallium experiments, we include the energy dependence of the cross section errors and assume 
full correlations among the cross section errors of the different neutrino sources that contribute at a specific energy. 
Specifically, we use the results from Ref. [ |56[ for the cross section errors. For continuum sources, we take the fractional 
error of the interaction cross section in the energy bin k to be 



AlnC Ga 



' max, mm 



(E k ) - a heat (E k )\ 



(Jbest(Ek) 



(A5) 



where cr max . min(£'fe) are the 3er upper and lower limit cross sections given in Tables III and IV of Ref. pq |. For the 
line sources, 7 Be and pep, we use the errors given in Eqs. (41) and (42) of Ref. pq| . To be complete, we have checked 
that the shapes of these neutrino lines (see Ref. |45|]) do not affect significantly the error estimates. The gallium 
contribution to the cross section error matrix is therefore given by 



dR 



dR 



fcl fc2 



d\nC{E kl ) dlnC{E k2 



-AlnC(H)AlnC(fc2) p klk2 , 



(A6) 



where p k \ k2 is the correlation matrix for the cross section errors of the different energy bins. For low energy neutrino 
sources, the dominant process is the transition to the germanium ground state and for higher energy neutrino sources 
(E > 2 MeV) transitions to excited states dominate. Therefore, we assume the cross section errors to be fully 
correlated between energy bins either below or above E = 2 MeV {p k \ k 2 — 1 for E k \,E k2 < 2 MeV or E k i,E k2 > 2 
MeV). We take the errors to be uncorrelated between one energy bin below 2 MeV and one energy bin above E = 2 
MeV (p klk2 = for E kl < 2 MeV and E k2 > 2 MeV, or E kx < 2 MeV and E k2 > 2 MeV). 

For chlorine, we take as fully correlated the cross section errors for the pp, pep, 7 Be, and CNO neutrinos. The cross 
section errors for the 8 B and hep neutrinos are uncorrelated with the errors for the lower energy neutrinos but are fully 
correlated with each other. We can neglect for chlorine the energy dependence of the chlorine to argon cross section 
errors because the uncertainty for the lower energy neutrinos is determined almost entirely from the ground-state 
to ground-state matrix element. (Forbidden corrections are unimportant for these low energy neutrinos.) For the 
8 B and hep neutrinos, the absorption cross sections are also dominated by a single (but different) transition, in this 
case the super-allowed transition [[[o| , and therefore we can also neglect the energy dependence for the higher energy 
neutrinos. We adopt the values of the averaged chlorine cross section errors given in Refs. To be explicit, the 

chlorine contribution to the cross section error matrix is given by 
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ct ckci ( c - s -) = £ J^iciifeci A In CftA In p iia (A7) 

il i2 

where pm2 = 1 for il, i2 = pep, CNO, 7 Be or il, i2 = 8 B, hep. Also, p ili2 = for il = pep, CNO, 7 Be and i2 = 8 B, hep. 
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